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Abstract 

We describe how to treat the interaction of travelhng electrons with locahsed vibrational 
modes in nanojunctions. We present a multichannel scattering technique which can be ap- 
plied to calculate the transport properties for realistic systems, and show how it is related to 
other methods that are useful in particular cases. We apply our technique to describe recent 
experiments on the conductance through molecular junctions. 



Introduction 

Electronic transport through nanoscale systems connected to external electrodes exhibits a num- 
ber of new features when compared to conduction through macroscopic systems. Primary among 
these is the importance of local interactions-including Coulomb interactions between the elec- 
trons, and scattering from localised atomic vibrations. Crudely speaking these effects are more 
important in nanoscale systems because the electronic probability density is concentrated in a 
small region of space; normal screening mechanisms are then ineffective. 

Perhaps the most extreme instance of nanoscale transport is where an electrical connection 
is formed through a single molecule. This paper describes how to approach electron transport 
in such systems theoretically. Our aim is to convince the reader that in many cases the local in- 
teraction is not a small perturbation which can be added in later — instead, it plays the principal 
role in determining the transport physics. In order to treat it, we need to understand the role 
of the local molecular properties — dominated by these local interactions — in scattering current- 
carrying electrons as they pass by. As well as determining the current flow, the interactions are 
then responsible for other important phenomena such as local heating. 

Such a situation is already familiar in the physics of transport through semiconductor quan- 
tum dots , where local Coulomb interactions produce Coulomb blockade (regions of very low 
conductance where changes in the charge state of the dot are energetically forbidden) and the 
Kondo effect (where the local spin of the dot electrons is screened by the conducting carriers) . 
But the situation in molecules, where the length scales are smaller and the electronic energy 
scales are correspondingly larger, is only now being studied. 

If the (three-dimensional) extent of an electronic wavepacket is L, the Coulomb interaction 
between electrons scales as ^ while the energy arising from short-range interactions with lat- 
tice vibrations scales as We might therefore expect the electron- lattice interactions to 
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dominate for single molecules where L is small, and indeed their signatures have been observed 
in several recent experiments with single molecules (or small numbers) El El El- These exper- 
iments include inelastic electron tunneling spectroscopy (lETS) of small molecules adsorbed on 
surfaces ^3,, studies of molecular-scale transistors @], and experiments on molecular junctions 
made of alkyl or 7r-conjugated molecular wires ^ or alkanedithiol self- assembled monolayers 
0. Such experiments show the possibility of exciting specific vibrational modes by injecting 
electrons; these excitations show up as sharp features in the differential conductance. 

It has been known for many years that coupling between electrons and molecular vibrations 
plays an important role in transport properties of bulk films of long conjugated molecules (such 
as conducting and semiconducting polymers); the transport then occurs by classical diffusion 
of polarons or solitons (HI • The new realisation is that such effects are also crucial in nanoscale 
transport through individual molecules, even although a classical diffusive model for the resulting 
composite excitations is no longer appropriate. This new realisation has coincided with an 
appreciation of the importance of local heating another manifestation of the coupling 

between electrons and phonons. 

Corresponding to this experimental progress, theorists have increasingly been studying elec- 
tron transport in the presence of local electron-phonon (e-ph) coupling in the molecule. In 
cases where the e-ph effects are weak (for example, in lETS of small molecules, where the elec- 
tron residence time on the molecule is very short), perturbation theory for the e-ph interaction 
can be used to interpret lETS spectra for model systems [5] and for more realistic systems 
including small molecules on surfaces |1(J| and molecular wires Going beyond perturbative 
approximations, the effects of e-ph coupling have been considered using scattering theory and 
Green's functions approaches |13l 1141 ITSl 1161 [T71 1181 using master and kinetic equations |2()| . 
and reduced density matrix approaches |21| . Another contribution including dissipation with 
a semiclassical treatment of the vibration was also developped More recently, approaches 
based on non equilibrium statistical physics have been developped for model systems I24j 
and realistic atomic and molecular wires 25j J6l |23 . 

In this paper, we describe the effects of the coupling between injected electrons and molecular 
vibrations (both treated on the quantum level) by using a multi-channel inelastic scattering 
technique. We show that it is important to treat such local interactions non-perturbatively, and 
how our approach relates to others. In the final part we study the effects of the temperature 
on the conductance properties and shed some light on the basic understanding of the features 
observed in lETS spectra. 

Note that in the following, we consider independent charge carriers interacting with sin- 
gle or multiple quantized vibrational modes. Calculations including the interaction between 
charges 28 or the interaction between charge and classical vibrations |29| have been considered 
elsewhere. 



Before presenting some numerical results, we describe our approach to inelastic transport and 
its relation to other approaches. 

Physical model 

Our model involves those delocalized electron states which simultaneously carry current through 
the molecule and interact with its eigenmodes of vibration. In the case of conjugated molecules 
these states are predominantly derived from the tt orbitals (though our method can deal equally 
with other types of state, and there are no restrictions on the geometry or dimensionality). 
If the molecular vibrations are approximately harmonic, the reference Hamiltonian Ho for the 
molecule is 
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where Hc\ is the purely electronic Hamiltonian and A labels the eigenmodes of vibration of the 
isolated molecule in its equilibrium (generally neutral) charge state; a\ creates {a\ annihilates) 
a quantum of energy uj\ in mode A. Note that Hd may, in principle, contain electron-electron 
interactions. 

Our fundamental approximation is to consider the scattering only of single electrons or 
holes from the molecule. We expect this assumption to be reasonable when the interval between 
carrier transmission events is much greater than the transit time through the molecule |17| . 
and provided that the applied bias does not cause significant fluctuations in the occupancy of 
the molecular eigenstates (i.e. provided one charge state of the molecule dominates during the 
transport process). 

For electron transport, if the equilibrium charge state contains N electrons and |0, N) is the 
A^-electron ground state, we construct an electronic basis \i,N + 1) = c||0, A^) where the {c|} 
create an electron in a sufficiently complete set of orbitals (for example, in all the low-lying 
TT-states of the molecule). For hole transport we would use \j,N — 1) = Cj\0,N). We now 
diagonalize Hd (keeping the phonon coordinates fixed, for the moment) within our restricted 
{N ± l)-electron basis to obtain a set of approximate {N ± l)-clcctron eigenstates. Using the 
electron creation (cj^) and annihilation (c„) operators and the energies e„, we can write 

Hci^^e„clcn- (2) 

n 

The approximate {N ± l)-electron eigenstates of the full Hq can thus be written as |n, {n\}), 
where n labels an electronic eigenstate and the {n\} are occupation numbers for the vibrational 
modes: 

\n,{nx})=ciY[{alr/V^. |0,iV,{0,}) (3) 

A 

and |0, N, {0\}) is the electronic and vibrational ground state for N electrons. 

The e-ph coupling term i/eph is taken to be linear in the phonon displacements and induces 
transitions between these electronic states: 

Heph^ ^ "fXnjnial + ax)clcm ■ (4) 
A,n,m 

The values of e„, uj\ and "fxnm are calculated from a suitable model of the isolated molecule. In 
our work on conjugated molecules ^7] we have used the Su-Schricffer-Heeger (SSH) jS] model 
within the harmonic limit. 

To obtain the transport properties, we perform a thought experiment in which we connect the 
left and right ends of the molecule to metallic leads and scatter a single incoming charge carrier 
(electron or hole). The coupling matrix elements are vl,r respectively; since we do not wish to 
focus on the properties of the leads we take the simplest possible model for them and consider 
one-dimensional semi- infinite chains (with on-site energy eL^R and hopping integrals I3l,r, giving 
dispersion relations e = el.r + 2f3L,RCOs{kL,R), where kL.R are dimensionless wavevectors) . We 
assume transport is purely elastic within the leads themselves; dissipation occurs in remote 
reservoirs, as in the standard Landauer picture for electron transport. The scattering states 
l^f) for a single incoming carrier are then expanded inside the molecule onto the eigenstates 
|n, {n>,}) of Hq. The single added carrier can be anywhere in the system and interacts with 
lattice vibrations only when inside the molecule. 

The transport problem is solved by mapping the many body problem onto a single-electron 
one with many scattering channels |13[ I15L [T7| ; each channel represents a process by which 
the electron might exchange energy with the vibration modes. For an initial mode distribution 
b = {m\} and an incoming electron from the left, the outgoing channels in the left and right leads 
are associated with energy-dependent reflection coefficients rah(e) and transmission coefficients 
iab(e), where a = {n\} is the final mode distribution. In the leads, the scattering states are 
propagating waves with amplitudes rat (refiection) and tat (transmission), and wave vectors 
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corresponding to the initial (ei„) and final (efin) electronic energies. By projecting out the leads, 
we work entirely in the molecular subspace to obtain the scattering state j^") by solving 

^G^H-i|vi/H) = |.,H). (5) 

Here uj is the conserved total energy; 

= <^in + ^mxoJx = '^tin + '^nxuJx- (6) 

A A 

The source term |sb(w)) is fixed by the incoming boundary conditions, while are the 

electron self-energies arising from the coupling of the molecule to the leads. The components 
of the scattering states then give the matrix elements of the retarded Green's function G""(tj). 
In practice, the linear system j^f) = G^\s) is solved for a finite-size molecular subspace, by 
truncating the occupation number above a maximum n™''^ in each mode. This is physically 
reasonable because the injected charge cannot populate infinitely many excitations. 

From the solution of Eq. jSj) , we can calculate any observable property, such as the expectation 
value of an operator A: (A) = I]{6} where W^'^iT) is the statistical weight 

of the distribution b = {toa} of modes at temperature T. The transmission probability Tab is 
given from the square of the transmission coefficients tab with the usual ratio of the electron 
velocities in the outgoing and incoming channels. In our model, tab o: {i — N\Gab\s) and the 
transmission probability Tab is written as [TU] 

? 2 

Tab{e&n,ein) = 4 -^sinfc^(ein) -r^sinfcf (efin) 
PL PR 

x\{^^N\G:,{u;)\^ = ^)\^ , (7) 

where {N\Gaij{uj)\l) is the matrix element of the Green's function G^ taken between the left 
side i — 1 and the right side i = A*" of the molecule and the vibration distributions before (6) 
and after (a) scattering. The factors w| ji/ Pl.r. sin k^'^ are related to the imaginary parts of 
the retarded self-energies E£ ^. 

Two concluding notes are in order. First, this procedure solves the problem non-perturbatively 
in the e-ph interaction. Second, we have described the method for the case where the channel 
structure of the leads is generated only by the vibrational excitations of the molecule, but it is 
straightforward to generalise it to the case of multiple spatial channels in the leads. 

Connection to other methods 

Here we etablish how the multichannel scattering technique is related to other methods based on 
two-particle Green's functions for non-interacting electrons and on non-equilibrium Green's 
functions |24l 1251 1261 ITfj . We will show below (for a simplified molecular model) that G^ can be 
reformulated to include the self-energies arising from the coupling of the molecule to the leads 
(^reads) ^'^'^ from the interaction of the electron with the vibrations {^eph)- 

Let us start with the generalisation of the Landauer formula for non-equilibrium interacting 
systems given in Ref.|32|. In the following we consider the single-site, single- vibrational mode 
(SSSM) model, which allows us to illustrate the essential physics without complicating the 
algebra. We first assume that the left and right leads are identical, and couple to the molecule 
via energy-independent hopping integrals. Then, the self-energies ^ (and their imaginary 
parts ^L,Fi) are proportional to each other. According to Ref . '32J , the current through the 
junction is then given by 

I^^Jduj ifLiuj) - fniuj)) ImTr{r(a;)G'-(c^)} , (8) 

where fL,R is the Fermi distribution of the left and right (non-interacting) lead respectively, 
F = Ti,Ti^/{Tl + Ffl) and C is the retarded Green's function of the molecule including the 
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self-energies arising from the coupling to the leads and from the interaction between particles. 
Here we consider non-interacting charge carriers coupled to vibrations of the molecule. 

Within the SSSM model there is only one energy level eo that couples to a single mode of 
frequency luq with a coupling constant 70. For this model, is the inverse of a tridiagonal 
matrix R in the subspace \xn) of the excitations n of the vibrational mode (n = is ground 
state of the mode) with diagonal elements Rn,n ^ to — eo ~ tiloq — T,^^^y{u}) and non-zero off- 
diagonal elements Rn,n+i ~ — 7ov?^ + l and Rn,n-i = — 7ov^- Here E^j^j(ijj) is the retarded 
self-energy of both leads for the channel containing n phonon excitations; it is calculated by 
taking into account the energy conservation condition Eq. , which reduces in the limit of very 
low temperature to uj — tin — ^fin + nujo. Furthermore, in this limit, only the matrix element 
of C between the phonon ground state Gqq = (xo|G''(a;)|xo) enters in the evaluation of the 
current. It can be expressed as a continued fraction: 



After some calculation, it can be shown that ImCpQ is actually related to the other non-diagonal 
elements {XnlG^lxo) of the Green's function and that finally one has 



Introducing this result in Eq.®, the current for the SSSM model is then expressed as the 
sum of the elastic and inelastic transmission probabilities Tab(e', e) obtained from Eq.Q. 

One can perform a similar derivation at a finite temperature T. The current is then de- 
termined from all the diagonal matrix elements Im (Xn|G''|xn) including the appropriate sta- 
tistical weight yVP^(r) for the n-th excitation of the vibration mode. Each matrix element 
Im {Xn\G^\Xn) is related to the sum of the elastic and inelastic transmission coefficients where 
the vibration mode is in the 71-th excited state before scattering, and in the m-th excited state 
afterwards. 

The derivations given above can be extended to general cases in which many electronic levels 
are coupled to many vibration modes. Then the R matrix which defines C is block-tridiagonal; 
the diagonal blocks correspond to purely electronic processes with constant phonon occupancy, 
while the off-diagonal blocks correspond to the e-ph coupling matrix "fxnm- The equivalent 
matrix element ({Oa}|G''|{Oa}) of Eq.® is then written as a matrix continued fraction in the 
corresponding basis. The same reasoning derived above holds for calculating the current in 
terms elastic and inelastic transmission probabilities Tab for all temperatures. 

Now, we show how to connect Ea. (|ll(l to the approach developed in For this, let us 

consider the wide-band limit for the leads, so the lead self-energies become independent of the 
energy and hence of the channel considered, and are purely imaginary Y^^^(lu) = iT = i{TL+Tf(). 
The imaginary part of Gqq then becomes Im Gqq{lu) = ^J2n \{Xn\G'^\xo)\'^ ■ It is convenient to 
work in this limit to circumvent the problems of band- width renormalisation when one introduces 
the transformation U that diagonalises H = Hq -\- Heph = ^qCqCq + ajQa^aQ + 70 (aj + ao)cJco. 
We introduce the basis set |xn) = U\xn), and obtain 



G[jo(cj) = - eo - S'|oj(w) - 
7o/(^ - eo - t^o - - 
2 7o/(cj - eo - 2cJo - ^\2}i^) - 
3 Jo/i^ - eo - 3wo - S{3}(^) - •••]"^ 



(9) 




(10) 




(11) 
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where eo = CQ — g'^ujo and g = jo/^a- Then, calculating the overlaps {xn\xi) between the original 
\Xn) and displaced harmonic states |xn) as in one ultimately finds 



°° 2n 

ImGSoM^re-2^ 

n=0 



X 



°° 2/ 

g 



j=o ^-^ ' 1=0 ^' ^ ^ ^0 ^ + + *r 



(13) 



Ea. ljl^ll is just the result derived in \IM\ for the transmission by starting from a two-particle 
Green's function description for a single resonant level coupled to a single vibration mode. 

It should be noticed that the equivalence between the two approaches is only valid in the wide- 
band limit for which the full polaron shift is obtained (i.e., eg shifted by the relaxation energy 
— 7q/wo). In other cases, the energy dependence of the self-energies S^^i plays an important 
role. Furthermore when the residence time of the electron in the molecule is not long enough 
for the vibration to respond fully to its presence, the full relaxation is not obtained and one 
has an intermediate polaron shift ^. However there are no such limitations in the multichannel 
scattering technique, which can treat the problem for the full parameter range (including strong 
e-ph coupling) — provided only that the fundamental assumption of single-carrier transport (see 
the discussion following Eq.(^) remains valid. 

Finally, in this section, we show how the multichannel scattering technique is related to 
more recent approaches based on non-equilibrium Green's functions 1^ 123 EHl EZ] • Such 
a non-equilibrium approach is more suitable when our assumption of a well-defined reference 
charge state breaks down. Let us first rewrite the retarded Green's function in Eq.Q as 

G^iuj) = - eo - S^o}(^) - Kpi^i^r' , (14) 

where the self-energy due to the interaction between electron and phonon is symbolically and 
recursively given by Egpj^(ci;) = G^{uj — ojq) ^. Note that by definition ^7], one has the 
following relation for the lead self-energy: S|pj(a; — nwo) = E|^j(ci;). 

In a many-body non-equilibrium Green's functions approach, the self-energies arising from 
the interaction between particles are obtained from a diagrammatic perturbation expansion 
of the interaction and by applying the rules for the time ordering and for the evaluation of 
products of double-timed operators on the Keldysh contour. In principle, for a coupled e-ph 
system, one should include processes to all orders in the interaction to calculate self-consistently 
the different electron Green's functions dressed by the phonons, as well as the phonon Green's 
functions dressed by the electrons. This is a tremendous task to achieve for realistic systems, and 
so far it has only been done almost exactly for model systems or within some approximations 
for atomic or molecular wires. Recent studies have been performed using the so-called self- 
consistent Born approximation |24l 1251 12()l l?7] , in which the retarded electron self-energy due 
to e-ph coupling is obtained from the sum of several contributions all written in the following 
form: Sgpj^_^-^(a;) oc / duj' {uj — lu')G'^ {lu') , where the superscripts X,Y = r,> represents 
the different types of Green's functions (retarded r, Keldysh greater >) for the electrons (G) 
and for the phonons (£>). 

In the following, we work with the undressed phonon Green's functions Df){Lj) for which the 
Keldysh greater component is Dq{lu) — ~i2TT{N{Lu)6{uj + ujq) + {N{uj) + l)6{u! — ujq)), N{u!) 



^It should be noted that no polaron shift is obtained if one makes a perturbative expansion of the e-ph coupling. 

^This is a closed form for a lowest order series expansion of the e-ph self-energy equivalent to the self-consistent 
Born approximation |24l 1251 1^127) . In fact the continued fraction expansion of E^pj^ includes a factor n at each level 
of the fraction as seen in Eq.@. Such factors arise from applying the creation (annihilation) phonon operator on the 
vibrational state \Xn-i{n)) in a multi-excitation process. The exact series expansion for EJIpi^ includes higher-order 
terms corresponding to multiple excitations, for example: 70 G^{iu — luo)G'^{w — 2LU())G^{iu — cuo). A similar result 
can be obtained by a linked cluster expansion approach. 
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being the Bose-Einstein distribution function. In the limit of low temperatures, simplifications 
arise because then N{u;) = 0. Furthermore by using the principle of causality and the analytic 
properties of C in the complex plane, it can be shown that the contribution Y^^^p^irri^) ~ 
—ijQ J duo' /2ti Dq{llj — u!')G^{uj') vanishes'^. The second contribution to the self-energy is given 
by / dui' /27r {ijj—lo')^' {ui') . At zero temperature, this is exactly the self-energy 7g G"'(w— 
Wo) derived previously in Eg. 114(1 . The final contribution to the self-energy is Egpj-j (w) cx 
/ duj' {Dq{uo')— Dq{uo = 0))G> {uj—u'). It involves the greater electron Green's function G> which 
provides information about the non-equilibirum density of unoccupied states of the molecule. 

is related to the corresponding self-energy S> = 2]jJ;^jg-)-S^pj^ and the retarded and advanced 
Green's functions G'''° via the kinetic equation G^ = G^Jy^G"" . The terms involving S^p^ are 
intrinsically of higher order in 0(79) than Y7^^T^{bj) — 7q G'^{uj — loq) and can be neglected for 
weak e-ph coupling. The other terms involving Sf^Jg^^^ are important; however their contributions 
become small when the electronic level coupled to the vibration is off-resonance with the Fermi 
levels of the leads (in the limit of zero to small applied bias). This is case for molecular wires 
having a substantial HOMO-LUMO gap, when the Fermi levels are pinned inside this gap. 

So, the multichannel scattering technique takes into account the most important contribu- 
tions to the e-ph self-energy. There are some terms arising from the Keldysh approach to trans- 
port that are not treated in the multichannel technique; however, as mentioned above, the gen- 
eral problem of the non-equilibrium transport is so complex that it has only been solved for model 
systems and using approximations for the Green's functions and self-energies [231 1241 1^ I26L I27| . 

Now that we have established how to relate the multichannel scattering technique to other 
widely used non-perturbative approaches for transport through e-ph coupled systems, we turn 
to some practical applications of the technique. 

Results 

In our previous studies on long conjugated molecules [181 117| . we have already shown that the 
longitudinal optic phonon modes are the most strongly coupled to the injected charges. We 
have elucidated the mechanisms of charge injection and transport in such molecular wires. The 
transport is associated with the formation and propagation of polarons in perfectly dimerized 
(semiconducting) molecular wires 17, . More complex mechanisms arise in the presence of mid- 
gap states involving the delocalisation of a soliton and polaron-soliton interactions il8) . 

Here we present new results for models of molecules in relation with recent experiments; 
namely we consider the effects of the temperature on the conductance peaks around the HOMO- 
LUMO gap for short molecular wires [SJ, and we point out the difficulties for understanding 
the features observed in lETS [51 E] in terms of single vibration mode analysis. In doing so, 
we provide a plausible explanation for the temperature dependence of the width and shape of 
conductance peaks from a model system. More importantly, we show that it is crucial to treat 
the coupling to all the relevant vibration modes simultaneously, instead of considering one mode 
at a time, in order to be able to identify the origin of the features in lETS spectra. 

Temperature effects on the conductance peaks 

Here we consider the experiments performed with a mechanically controlled breakjunction at 
room temperature and at low temperature (T « 30K) for polyphenylene-based molecular wires 
|34|. The conductance curves show a peak just above a region of very low conductance at small 
bias (which might be a HOMO-LUMO gap or Coulomb gap). At low temperature, the peak 
is asymmetric with a maximum width of ~125 meV. At room temperature, this peak has a 
reduced amplitude and is more symmetric with a bell-like shape of maximum width 300 meV. 
It is clear that such a temperature dependence cannot be explained by the broadening of the 

^This is always true for electron-hole-symmetric systems. It is also true when the system is not too far from 
equilibrium. 
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T=0 room T 




CO [eV] CO [eV] 

Figure 1: Total effective transmission Tgfj versus energy u for the SSSM model (eo ~ 0.5 eV, loq 
80 meV, 70 ~ ujq/2). Left panel: Very low temperature limit: the transmission peak (red line) 
is asymmetric towards higher energy because of the presence of phonon side-bands corresponding 
to phonon emission. The resonance peak obtained in the absence of e-ph coupling (green line) is 
symmetric around eo- Right panel: At room temperature, both phonon emission and adsorption 
processes are available. The transmission peak (red line) is more symmetric around cq and its 
amplitude is reduced. In both panels, the blue lines represent the transmission obtained for a 
strongly reduced broadening in order to show explicitly the phonon side-bands. 

Fermi sea of the electrodes: a broadening over kT is not enough to explain the widening of 
the conductance peak and its change of shape. We therefore assume that there are fluctuations 
within the junction or at the molecule-electrode contacts, and that these fluctuations are related 
to low frequency vibration modes. It is difficult to identify exactly the nature of such modes, 
so we consider for simplicity a SSSM model in which the single level eg is coupled to a low 
frequency vibration mode {ujq « 80 meV) via a coupling constant 70 « wo/2. 

We perform calculations for this model according to the prescriptions given in the previous 
section at very low temperature and at room temperature. The total effective transmission 
Teff (w) is obtained as the sum of all the elastic and inelastic contributions to the transmission 
and is shown in Fig. ^ for the zero temperature limit {kT = 0) and for room temperature 
{kT 0^0/2). In the absence of e-ph coupling, the transmission is given by a single symmetric 
resonance located around eo, as expected. In the presence of e-ph coupling, the shape of the 
resonance is modified due to the presence of phonon side-bands. For very low temperatures, the 
transmission peak is asymmetric towards higher energy. The phonon side-band peaks are located 
at energies above eq because only phonon emission processes are allowed at low temperatures. 
At room temperature, both phonon emission and absorption processes become possible and the 
phonon band peaks appear for energies above and below the molecular level.'' The transmission 
peak is then more symmetric around eq. Futhermore, its amplitude is reduced compared to the 
peak at T = 0, as observed experimentally. 

^This lineshape appears similar to that obtained for very low temperature and very strong coupling |33|. but the 
physics is quite different because absorption processes are important. 
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Figure 2: Total effective transmission Tgg versus energy u) for a molecular wire of length A'' = 155. 
The transmission is calculated by including a single vibration modes one at a time. Optic mode 
= 112.1 meV (red line) showing the most important polaron shift, tox^ = 129.9 meV (green 
line), = 136.8 meV (blue line). The transmission in the absence of e-ph coupling is also shown 
(black line) . The vertical bars represent (for the same color) the eigenvalues of the isolated molecule 
Hamiltonian including the coupling of the single mode A to the many molecular levels. 

It is interesting that such a simple model provides results in qualitative agreement with 
the conductance measurements performed in Ref . . Such a mechanism therefore provides a 
plausible explanation for the temperature dependence of the conductance peaks in molecular 
junctions, although we cannot rule out other possible explanations such as those given in Refs. 

Multi-vibration mode effects 

We now consider a more realistic molecular wire built from an odd number N of monomers 
as obtained from the SSH model. In this case, there is a mid-gap state in the HOMO-LUMO 
gap. In the following, we couple all the unoccupied tt states to the low energy optic modes as in 
Eq. ^ . The actual spatial and frequency structure of these modes is very important in allowing 
charge to propagate through the molecule — in contrast to previous studies on simpler model 
systems [T5) . 

As a generalisation of the SSSM model, there will also be phonon side-bands in the transmis- 
sion. The phonon side-band peaks related to the mid-gap state will appear inside the HOMO- 
LUMO gap. Such peaks also exist in the absence of the mid-gap state, i.e. for wires of even 
length N; however, there they are hidden in the background of the resonances from the molecular 
conduction band. 

The transmission curves for a wire of length N = 155 in which the molecular levels interact 
with a single vibration mode one at a time are shown in FigI21 Several different optic modes 
(wai,A2,A3 = 112.1, 129.9, 136.8 meV) are considered. The total effective transmission Tcff shows 
phonon side-band peaks in the HOMO-LUMO gap, separated as expected by the corresponding 
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Figure 3: Total effective transmission T^s versus energy u) for a molecular wire of length N = 155. 
The transmission (magenta line) is calculated by including the coupling to the three optic modes 
shown in Fig|21 New spectroscopic features appear in the transmission. The vertical bars give 
the eigenvalues of the isolated molecular Hamiltonian with many levels coupled to the three optic 
modes. The transmission in the absence of e-ph coupling is also shown (black line) as well as the 
transmission obtained for a coupling to a single mode (thin solid line) as in Fig. |2l 

energy u!\. Up to small energy shifts due to real part of the self-energies, the peaks in the 
transmission correspond to the eigenvalues of the coupled e-ph Hamiltonian of the isolated 
molecule, as shown by the vertical bars in Fig[21 For some eigenvalues, there are no resonances 
in the transmission. There are several possible reasons for this: (i) the phonon side-band 
peaks associated to a given electronic transition have an exponentially decaying weight for 
large numbers of excitation quanta, so some of these peaks may disappear in the tail of other 
resonances; (ii) there might be destructive quantum interferences for the corresponding process^; 
(iii) the corresponding eigenstates are such that the matrix elements (^|G^j|l) in Eq.Q are 
small or vanishing. 

Note that such calculations also reveal which optic phonon mode contributes the most to 
the polaron shift (shift towards lower energy of the first resonance above the gap due to the 
e-ph coupling). In the present case, the optic mode of energy uj\-^ = 112.1 meV gives the most 
important polaron shift. 

However for realistic molecules, the e-ph coupling matrix is non-diagonal and at the same 
time the injected charge may couple simultaneously to different vibration modes. Calculations 
including simultaneous coupling to the three optic modes mentioned above were performed. 
New spectroscopic information is obtained, as shown in Fig. 01 The transmission is not simply 
a superposition of the contribution of the three different modes shown in Fig. |21 instead, new 
peaks appear. These features can be interpreted in terms of quantum interferences between 
the different elastic and inelastic electron paths as well as in terms of beating effects. All 

''Such effects appear more often when several electronic levels are coupled simultaneously to several vibration 
modes. 
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these features in the transmission correspond to the eigenvalues of the coupled e-ph molecular 
Hamiltonian H, provided coupling of the many molecular levels to the three different vibration 
modes is included. A detailed analysis of the origins of such features is beyond the scope of the 
present paper; however, these effects show that the general form for the e-ph coupling Eq.l^J in 
realistic systems creates features in the transmission that cannot be obtained in terms of single 
vibration mode analysis. 

This is a very important result when one considers the interpretation of recent lETS exper- 
iments in molecular junctions^. Our results shown that special care is needed when assigning 
features in the transmission or the current to specific vibration modes of the molecular junc- 
tion. Calculations for realistic three-dimensional geometries in molecular junctions are currently 
under investigation. 
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